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Abstract 

The effect of the finite height of a cylindrical lead salt nanorod on its electronic structure is 
studied within the effective mass formalism. It is demonstrated that for most practical purposes it 
can be accounted for by sampling the subband dispersion dependencies of the infinite cylindrical 
quantum wire of the same radius at wave numbers k z = im z /H, where H is the nanorod height 
and n z is an integer. However, under certain conditions the otherwise degenerate electron energy 
levels will repel. A detailed account of these conditions is provided. 
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In recent years a significant progress has been made in solution-based synthesis of semi- 
conductor nanostructures. As a result, a new class of nanostructures called nanorods has 
emerged jl|. The nanorods can be considered as an intermediate case between the quasi- 
zero-dimensional (OD) quantum dots and quasi-one-dimensional (ID) quantum wires and 
allow one to investigate the variation in material properties in the transition from OD to 
ID. They can be thought of as cylindrical semiconductor structures having diameter on the 
order of several nanometers and height up to several tens of nanometers. 

Recently nanorods of lead salt semiconductors have been synthesized by a number of 
research groups [2j-|5]. Lead salt semiconductor compounds (PbS, PbSe, PbTe) are charac- 
terized by a large bulk exciton Bohr radius, in the range of 20 to 46 nm. Thus, in lead salt 
nanorods of a relatively short height, the strong confinement regime of charge carriers can 
be achieved in all three dimensions. 

In Ref. [6] we have studied electronic structure of cylindrical lead salt quantum wires 
within the effective mass approximation. We found that the electron energy subbands 
of the conduction and valence bands in these quantum wires have monotonous dispersion 
dependencies. At first glance, in this situation the effect of finite dimension of the structure 
along the nanorod growth direction on the electron energy structure could be accounted for 
by substitution of the discrete values of wave numbers, im z /H (where H is the nanorod 
height and n z is an integer) into energy dispersion dependencies for various subbands of 
the quantum wire. It turns out, however, that such procedure does not account for some 
qualitative changes in electron energy structure caused by the finite height of the nanorod. 
In this paper we will first describe these changes using a numerical solution of the effective 
mass equations, and then explain them based on the analytical dispersion equation derived 
in Ref. 6| for the case of a quantum wire. 

The conduction and valence band extrema in lead salt semiconductors occur at the L- 
points of the Brillouin zone leading to the high valley degeneracy. Bartnik et al. have studied 
solutions of the effective mass equations customized for each of the four L- valleys for different 
quantum wire growth directions [7j. They showed that the valleys remain degenerate, within 
the effective mass approximation, when the growth occurs along the (100) direction. It is not 



infrequent that lead salt nanorods grow along this crystallographic axis 



5|, lTj] . The valley 



2 



effective masses in this case still remain anisotropic. In what follows we will neglect this 
anisotropy and assume that the nanorods grow along the (100) direction. Then the electron 
spectrum near the L-point taking into account only the two closely lying conduction and 
valence bands can be described by the spherical Dimmock model [6 
electron wave function is written as 



In this model the 
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where |L 6 ) and |Lg~) describe the Bloch functions while u(r) and v(r) are the spinors slowly 
varying with coordinates and satisfying the equation 
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Here ap ((3 = x,y,z) are the Pauli matrices, a c , a v , E g , and P are parameters of the 
model and E is the electron energy. Eq. ([3]) is written using the system of units where 
h = ttlq = e = 1 (mo and e are the mass and the charge of a free electron). Inclusion of the 
valley anisotropy would amount to introduction of longitudinal and transverse counterparts 
for the parameters P, a c , and a v Q]. 

In order to solve Eq. (T5]) numerically for boundary conditions of cylindrical symmetry, 
one can first write down the Hamiltonian ([3]) in the basis of four cylindrically symmetric 
bispinors of the form 
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Here R is the radius of the cylindrical nanostructure, xi^ n is the nth zero of the Bessel 
function, Jl(x), and 



V2 



is the normalization coefficient. All these bispinors are characterized by the quantum number 
M, the projection of the total angular momentum onto the growth direction, z, and vanish 
on the cylindrical surface of the nanostructure. 

For a cylindrical quantum wire the function f(z) should be chosen in the form of a plane 

wave 

/(*) 
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If the index n runs from 1 to n max , then the Hamiltonian fl3]) written in the chosen basis, for 
given quantum numbers M and k z , represents a 4n max x ^n max matrix. If one is interested 
in a certain number of subbands characterized by the given value of the total angular mo- 
mentum projection, M, it is enough to set n max equal to that number. One can check that 
numerical diagonalization of the resulting Hamiltonian leads to electron energy dispersion 
dependence coinciding (with a very high precision) with that obtained from solution of the 
dispersion equation derived analytically and given by Eq. (25) of Ref. 6]. The numerical 
procedure is very robust and for practical calculations it would be our first choice. 
For a nanorod of height H the function f(z) can be chosen in the form 



, / 2 . im z z 



(4) 



If the index n z runs from 1 to n™ ax then the Hamiltonian ([3]) , for a given quantum number 



M, represents in the chosen basis a An max n™ ax x An 



max '^z 



matrix. We have numerically 
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FIG. 1: (Color online) Energies of levels of size quantization described by the quantum number 
M = 1/2 in a PbSe nanorod of radius R = 20 A and height H as functions of H/R (solid lines). 
Dashed lines: energy dispersion Ei/2 >ric (k z ) for subbands of a PbSe cylindrical quantum wire of 
radius R = 20 A for quantum numbers M = 1/2, n c = 1,2 [6[ calculated at k z = 7m z /H. 

diagonalized this matrix for nanorod radius R = 20 A and nanorod heights in the range 
from 10 A to 200 A. We used the same material parameters as in Ref. [6| (corresponding 
to PbSe nanorods) and the values of n max = 3, n™ ax = 30. In Fig. 1 by solid lines are 
shown the dependencies of the first 10 eigenvalues occuring within the conduction band and 
corresponding to the quantum number M = 1/2 on the height to radius ratio, H/R. By 
dashed lines in Fig. 1 are shown the quantum wire subband energies |6[ E 1 / 2 ,n c (k z ) calculated 
at k z = irn z /H for a PbSe quantum wire of radius R = 20 A. 

The fact that the solid lines in Fig. 1 closely follow the dashed ones suggests that the 
quantum number n z and the wave function (jlj) represent a good zero-order approximation 
to describe the motion along the nanorod growth direction, z. On the other hand, one can 
see that the dependencies of nanorod energy levels presented in Fig. 1 demonstrate both 
crossings and anti-crossings (or avoided crossings). In order to reveal the underlying pattern 
some regions of Fig. 1 are presented in Figs. 2 and 3 on an enlarged scale. These figures 
reveal that the level crossings and anti-crossings form a checker-board pattern, although the 
anti-crossings become less pronounced as the nanorod height increases. The same situation 
occurs for other values of M, as demonstrated in Fig. 4 for M = 3/2. In Fig. 5 the levels 




FIG. 2: (Color online) Same as Fig. 1 but on an enlarged scale. 
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FIG. 3: (Color online) Same as Fig. 1 but on an enlarged scale. 



for \M\ = 1/2 and \M\ = 3/2 are plotted together to account for lowest conduction band 
energy levels in nanorods. 

The occurrence of the checker-board pattern of crossings and anti-crossings in the height 
dependence of electron energy levels in lead salts nanorods can be understood if one analyzes 
the equation describing electron energy dispersion in lead salts quantum wires 6|. This 
equation is given by Eq. (25) of Ref . 6] . For a given value of the quantum number M it has 
the following structure 

k 2 z A + BC = 0. (5) 
When k z = Eq. (JSJ) breaks into two separate equations B = and C = describing 
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FIG. 4: (Color online) Same as Fig. 1 but for M = 3/2. 




FIG. 5: (Color online) Energies of levels of size quantization for the quantum number \M\ =1/2 
(solid lines) and \M\ = 3/2 (dashed lines) in a PbSe nanorod of radius R = 20 A and height H as 
functions of H/R. 

quantum states of opposite parities at the apexes of ID subbands 6]. When k z ^ the two 
solutions mix. This mixing is fairly weak and can be considered as a perturbation of the 
subband dispersion relations described by the equations B = and C = 0. The structure 
of Eq. ([5]) suggests that this perturbation is linear in k z . In the space of envelope wave 
functions k z can be considered as the operator of the z component of the linear momentum 
in the momentum representation. The same operator in the coordinate representation is 
proportional to the first derivative with respect to z, or d/dz. The functions (jlj) are either 



odd or even with respect to the transformation z — >■ H — z. The operator d/dz can only 
have non-vanishing matrix elements between the functions of different parity (under the 
transformation z — > H — z). Thus, the levels originating from the lower quantum wire 
subband characterized in the zeroth order by odd (even) n z cross with the levels originating 
from the upper quantum wire subband characterized by odd (even) n z and anti-cross with 
the levels characterized by even (odd) n z , provided that these subbands have the same 
quantum number \M\. Note, however, that, if the states in the first (n c = 1) quantum wire 
subband with a given \M\ are characterized at k z = by the equation B = so will the 
states from the third (n c = 3) subband. Therefore, the electron energy levels in nanorods 
originating from the first and the third quantum wire subbands characterized by the same 
quantum number \M\ will always cross, and the checker board pattern will be broken for 
higher energies. This observation has been verified by numerical calculations. 

To summarize, we have demonstrated that for most practical purposes the effect of the 
finite height of a lead salt nanorod on its electronic structure can be accounted for by 
sampling the subband dispersion dependencies of the infinite cylindrical quantum wire of 
the same radius at wave numbers k z = 7rn z /H, where H is the nanorod height and n z is an 
integer. This approximation can fail if, for a given nanorod height, some of the energy levels 
originating from different quantum wire subbands characterized by the same absolute value 
\M\ of the total angular momentum projection onto the quantum wire axis are close to one 
another. The levels would repel if the corresponding states are characterized by the numbers 
n z of different parities and originate from subbands characterized by the same value of \M\ 
and main quantum numbers of opposite parities, giving rise to the levels splitting. 
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